Wannier-based definition of layer polarizations in perovskite superlattices 
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In insulators, the method of Marzari and Vanderbilt [Phys. Rev. B 56, 12847 (1997)] can be 
used to generate maximally localized Wannier functions whose centers are related to the electronic 
polarization. In the case of layered insulators, this approach can be adapted to provide a natural 
definition of the local polarization associated with each layer, based on the locations of the nuclear 
charges and one-dimensional Wannier centers comprising each layer. Here, we use this approach to 
compute and analyze layer polarizations of ferroelectric perovskite superlattices, including changes 
in layer polarizations induced by sublattice displacements (i.e., layer-decomposed Born efi'ective 
charges) and local symmetry breaking at the interfaces. The method provides a powerful tool for 
analyzing the polarization-related properties of complex layered oxide systems. 
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PACS numbers: 77.22.-d, 77.22. Ej, 77.80.-e, 77.84.Bw 

Multicomponent superlattices based on the ABO3 per- 
ovskite structure have received much attention recently 
due to the exciting properties they possess as multifunc- 
tional materials (see Ref. Q and references therein) . Ex- 
perimental studies using modern layer-by-layer epitaxial 
growth techniques have gone hand in hand with accurate 
first-principles calculations that have helped to interpret 
experimental results and to guide the search for super- 
lattice compounds with tailored properties. For exam- 
ple, a compositional perturbation that breaks inversion 
symmetry was predicted 2] to allow tuning of the dielec- 
tric and piezoelectric response, as confirmed later when 
such superlattices could be grown experimentally 0, . 
In addition, strained-layer superlattices can show a sub- 
stantial enhancement of spontaneous polarization; such 
effects have been observed and analyzed using first- 
principles calculations ^ . 

One issue that has received much attention theoreti- 
cally is how to quantify the concept of local polarization. 
This can be very useful in isolating the contributions of 
constituent layers to dielectric and piezoelectric proper- 
ties 0, 111 , separating the effects of factors such as epitax- 
ial strain and applied electric fields 0,0i0,0|' ^^'^ under- 
standing the enhancement or suppression of spontaneous 
polarization 0, Q- A local description is also essential 
for characterizing and understanding interface contribu- 
tions to such properties. Among previously proposed lo- 
cal approaches H, , that of Meyer and Vanderbilt 
has been one of the most commonly used 0, S ■ Based 
on a linear approximation involving effective charges and 
small ionic distortions from a higher-symmetry nonpolar 
reference structure, this simple model captures the essen- 
tial physics and provides a semi-quantitative description 
useful for understanding many aspects of the behavior 
of multicomponent superlattices. However, as we shall 
discuss, it is neither exact nor unique. 

A first-principles method for identifying local dipolcs 
and computing their dipole moments in an extended sys- 
tem has been proposed in Refs. 0, Q|, based on ex- 



pressing the electric polarization in terms of the centers 
of charge of Wannier functions (WF) [l^ that are 
maximally localized along the direction of interest [l^ . 
The method was successfully applied to analyze the per- 
mittivity of ultra-thin Si-Si02 heterostructures |lfllllj. 

Here, we further develop a closely related method 
and apply them to ferroelectric perovskite superlat- 
tices, where the polarization, particularly that along the 
growth direction z, is of central physical importance. We 
introduce a WF-based expression for the "layer polariza- 
tion" along z associated with each charge-neutral AO or 
BO2 layer in an (001) superlattice built from II- VI ABO3 
perovskites such as BaTiOa, SrTiOa, and PbTiOa. Un- 
like the approach of Ref.|9|, the present one is exact (i.e., 
the sum of LPs relates exactly to the total supercell po- 
larization) and is entirely free of arbitrary choices in its 
implementation. We will present examples showing how 
this approach naturally provides an insightful local de- 
scription of the polarization behavior of perovskite su- 
perlattices, both at zero electric field and under nonzero 
electrical bias, and in particular yields valuable informa- 
tion about the highly localized atomic and electronic re- 
arrangements at the interfaces. 

The modern theory of polarization jl2j is routinely 
used to compute the polarization of a crystal as a sum of 
ionic and electronic (Berry phase) contributions. In the 
Wannier representation, this takes the form 
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where r and m run over ion cores (of charge Qr located 
at Rt) and Wannier centers (of charge — 2e located at 
fm), respectively, in the unit cell of volume V. In the 
case of a II- VI perovskite superlattice, one may hope to 
decompose the system into neutral "layers" (that is, AO 
or BO2 subunits) and define a "layer polarization" (LP) 
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in which the sums are restricted to entities belonging to 
layer j. Here S is the basal cell area and we are now 
focusing only on z components. The LP pj thus defined 
has units of dipole moment per unit area, and the total 
polarization, with units of dipole moment per volume, is 
exactly related to the sum of LPs via Pz = c^^J2jPj 
where c — V/S is the supercell lattice constant along 
z. For such a decomposition to be meaningful, we need 
(i) to resolve the arbitrariness associated with the posi- 
tions of the Wannier centers, and (ii) to be satisfied that 
the Wannier centers can be assigned to layers without 
ambiguity. We shall show below by example that (ii) is 
satisfied for the systems of interest, and thus we next 
turn our attention to issue (i). 

As is well known, the Wannier centers f,„ = 
(Wm|r|Wm) are not unique because the electronic struc- 
ture is invariant to unitary rotations among the WFs 
(corresponding, e.g., to different choices of phases of the 
Bloch functions in k-space) . Marzari and Vanderbilt 0| 
introduced a method for obtaining a unique set of WFs 
by choosing the ones that minimize the sum of second- 
moment spreads (spatial variances) of the WFs. In a 
three-dimensional system, this involves finding a best 
possible compromise between minimal spread inx,y, and 
z directions, and an iterative procedure is needed to find 
this compromise solution. 

Here, we are interested only in polarizations along z, 
and can limit ourselves to minimizing the spread only 
in that direction [lol Moreover, we can use a hy- 
brid representation of the electronic ground state that 
is Bloch-like in x and y and Wannier-like only along 
z We start from a conventional band-structure 

calculation carried out on a mesh of reciprocal points 
k = (kx,ky,kz) and adopt the relabeling g = {kx,ky) 
and k — kz- That is, each 2D vector g labels a string of 
J k-points running along the z direction with separation 
6 = 27r/ Jc. Our task is then to transform the Bloch func- 
tions IV'g.n/c) (n = 1, . . . , iV) into hybrid WFs ^jj | Wg,„i) 
(to = 1, . . . , N) via a ID Wannier transform, where is 
the number of occupied bands. We will then let z^ in 
Eq. (0) be the average of Zm(g) = (Wg,m|z| Wg^™) over 
the 2D mesh of g points. Since the ID Wannier trans- 
form is done independently at each g, we drop the g label 
in the following paragraph. 

The case of maximally localized WFs in ID was treated 
expHcitly in Sec. IV.C.l of Ref.l?. There, it was shown 
that the WFs that minimize the spread functional are 
identical to the eigenfunctions of the projected posi- 
tion operator PzP, where P = X^nfc is the 
band projection operator. It was also shown how they 
could be obtained from a parallel-transport based con- 
struction using the singular value decomposition (SVD) 
of the overlap matrices between neighboring k-points, 
Mmi = {umk\un,k+b) , whcrc Umk IS the periodic part 
of the Bloch function ipmk- 

The SVD is M = VTW'^ where V and W are unitary 
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FIG. 1: (Color online.) Dispersion of WF center positions 
for a IBT/IST superlattice as a function of g = {kx, ky). 



and E is positive real diagonal. For small 6, E approaches 
the unit matrix, and omitting E to write M = UV^ can 
be regarded as a way of constructing a purified version 
of M that is exactly unitary. Then A = HLi ^^''^ 
defines a global unitary matrix describing the parallel 
transport of the states on the k-point string, and its uni- 
modular eigenvalues Am define the Wannier centers via 
Zjn — (— c/27r) Imln A„i. Note that no iterative proce- 
dure is required; these ID Wannier locations can be ob- 
tained by a straightforward small-matrix diagonalization. 
A procedure that is similar in spirit, but slightly different 
in detail, has recently been proposed elsewhere 

In Fig. we present the resulting values of Zm(g) for 
an ab-initio calculation on a 10-atom tetragonal super- 
cell composed of alternating SrTiOs (ST) and BaTiOs 
(BT) units, which we refer to as a IST/IBT superlat- 
tice. We did all calculations using the ABINIT code 
'w], which implements density- functional theory within 
the local-density approximation (LDA) .17.|. We adopted 
Teter norm-conserving pseudopotentials [l8| for which 
the valence states are (5s5p6s) for Ba, (4s4p5s) for Sr, 
(3s3p3c?4s) for Ti, and (2s 2p) for O. We used a plane- 
wave energy cutoff of 45 Ha, a 6 x 6 x 3 Monkhorst-Pack 
self-consistency mesh, and a 12 x 12 x 3 bandstructure 
mesh. We assumed perfect epitaxial growth of the su- 
perlattices on a cubic ST substrate having a theoretical 
equilibrium lattice constant of 7.265 bohr and tetragonal 
P4mm symmetry. 

The key feature visible in the ^(g) dispersion relation 
in Fig.^is that the WF centers separate quite naturally 
into distinct layers as anticipated. The ID Wannier po- 
sitions z are almost independent of g = {kx,ky), and 
there are robust gaps between layers. Moreover, we find 
eight Wannier centers in each BaO or SrO layer and 12 
in each Ti02 layer (four for each cation semicore shell 
and four for each oxygen 2s2p shell), so that the layers 
are neutral as expected. All this demonstrates that the 
proposed Wannier-based approach does indeed lead to a 
natural and robust decomposition into easily identified 
neutral layers. 
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FIG. 2: (Color online.) Local polarization profile of 1ST/2BT 
supercell, from effective charge approximation based on A 
centered (open square) or B centered (open diamonds) anal- 
ysis, and from layer polarization analysis (filled circles). The 
overall supercell polarization is 0.22 C/m^. 

It is then straightforward to define the LP pj associated 
with each layer according to Eq. Q . For comparison with 
other definitions, we also introduce the corresponding lo- 
cal polarization Pj — Pj/cj having the correct units of 
polarization (dipole per unit volume), where Cj is chosen 
as half the distance between the two neighboring cations. 
In Fig. 121 we show the local polarization profile calcu- 
lated in this way for the case of a 1ST/2BT superlattice 
(15-atom supercell). The details are the same as for the 
calculation of the IST/IBT superlattice, except that we 
use a 6 X 6 X 2 k-point sampling. We compare our re- 
sults with the ones obtained from the commonly applied 
approximate scheme 9j in which the local polarization 
is estimated by multiplying the Born effective charges of 
the atoms in a unit cell layer by their displacements rel- 
ative to a reference structure. The effective charges are 
obtained from linear response calculations in the ferro- 
electric ground state. By its very nature, this approx- 
imate scheme |^ has only half the spatial resolution of 
our new scheme, since it applies only to entire ABO3 
cells. (Note that it cannot easily be extended down to 
the resolution of AO and BO2 layers because the sum of 
Z* values in such a layer does not vanish, so that the 
definition would depend on choice of reference structure 
in an unsatisfactory way.) 

The results shown in Fig. |21 are consistent with the 
findings of previous theoretical studies Q showing that 
the SrTiOs portion of the supercell becomes polarized to 
almost the same degree as the BaTiOa portion. How- 
ever, the improved resolution associated with the new 
approach is also clearly evident in the figure. For exam- 
ple, one can now see that the polarization tends to be 
larger in the AO layers than in the Ti02 layers (see next 
paragraph). Moreover, our new approach is free of three 
limitations of the approximate one 0, IE 0- First, 
we avoid the choice between an A- or B-centered anal- 
ysis. Second, we do not have the problem of choosing 
an arbitrary local reference structure as the basis for the 
definition of the atomic displacements. Third, our LPs 
Pj sum to give the exact total polarization of the entire 



supercell, whereas the approximate ones do not. 

The Born effective charges Z* , defined as the first 
derivatives of polarization with respect to atomic dis- 
placements, describe the dynamics of the charge transfer 
induced by such displacements. We now illustrate how 
the LP concept can be used to decompose the Z* for 
an atom in one layer into contributions from neighboring 
layers 19]. This is demonstrated in Tabled for the case 
of a supercell of tetragonal bulk BT that has been tripled 
along [001] (3BT supercell). Each of the four symmetry- 
inequivalent atoms (Ti and 0|| in a Ti02 layer and Ba 
and O^ in a BaO layer) was displaced along z in turn, 
and the changes in all six LPs in the supercell were com- 
puted. For Ba, Oj^ and 0||, the induced polarizations 
are dominated by contributions from the same atomic 
layer, at the level of around 45%. In contrast, for the Ti 
atom, the contributions from the first neighboring lay- 
ers are almost as large as from the layer itself. This is 
consistent with the well-known role of the Ti(3c?)-0(2p) 
hybridization in giving rise to the anomalous Born effec- 
tive charges in these perovskites ^J, and the fact that 
the WFs that embody this hybridization reside on the O 
atoms. Thus, a motion of the Ti atoms along [001] modu- 
lates this hybridization and shifts the centers of the WFs 
residing on the neighboring BaO layers. This effect also 
helps explain why the LPs for AO layers are larger than 
for Ti02 layers in Fig. [3 

We further illustrate the utility of the LP analysis by 
considering its behavior in a macroscopic electric field 
£ :2li applied along [001]. We used a constrained- 
polarization mapping technique |22| generalized to in- 
clude volume relaxation 23] to find the minimum-energy 
configuration for each given polarization. The resulting 
LPs Pj vs. total polarization are shown for the 1ST/2BT 
supercell in the left-hand panels of Fig. O We see that 
each LP is roughly linear in the total polarization for all 
six layers, but with nonlinearities appearing at large val- 
ues of polarization. For each data point, we extracted the 
corresponding macroscopic electric field value, and plot- 
ted the LPs against this field in the right-hand panels 
of Fig. 131 The results show a strongly nonlinear depen- 



TABLE L Layer decomposition of the [001] Born effective 
charges in a 3BT supercell. Total effective charges are given 
in the last row. 





Ti (IB) B 


a (lA) 0| 


(lA) 0_ 


L (IB) 


BaO (lA) 


1.433 


1.268 - 


-2.448 - 


-0.225 


Ti02 (IB) 


1.872 


0.148 - 


-0.231 - 


-0.930 


BaO (2A) 


1.262 


0.434 - 


-1.027 - 


-0.191 


Ti02 (2B) 


0.619 


0.296 - 


-0.542 - 


-0.216 


BaO (3A) 


1.211 


0.435 - 


-1.046 - 


-0.348 


Ti02 (3B) 


0.636 


0.191 - 


-0.264 - 


-0.217 


Z* 


7.033 


2.772 - 


-5.557 - 


-2.127 
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FIG. 3: (Color online.) Layer polarization as a function of 
(a) total polarization of the supercell, and (b) macroscopic 
electric field in the supercell, for six consecutive layers in the 
1ST/2BT supercell. Labeling of layers follows Fig. 1. 



dence, typical of that found for total polarization as a 
function of electric field in ferroelectric materials. 

The most striking features seen in Fig. O are for the 
two interface layers (Ti02 layers IB and 3B). They show 
local breaking of inversion symmetry; that is, the LP as 
a function of macroscopic field does not pass through the 
origin, at which for the system as a whole the symmetry 
of the PA/mmm space group requires P = and 5 = 0. 
One of these interface Ti02 layers has a nearest-neighbor 
SrO layer above and a BaO layer below, and the other 
vice versa, the two interface layers being related by a mir- 
ror symmetry. The LP approach give us much more pre- 
cise information about the response of these interfaces to 
applied fields than could be obtained from an analysis of 
either the total polarization, or of the local polarizations 
as previously defined (involving a smearing over three se- 
quential atomic layers). We expect that this method of 
analysis of interface layers will be invaluable for identi- 
fying the interface contributions to the properties of su- 
perlattices with more than two components, particularly 
those with globally broken inversion symmetry ^ .231] . 
Similar considerations apply to BaO layers 2A and 3A, 
which also see an environment of broken inversion sym- 
metry. For these, however, the symmetry breaking enters 
only at the level of second- neighbor layers, so the effects 
are smaller in magnitude. 

In summary, we have introduced a definition of the 
layer polarization (LP) in a multicomponent perovskite 
superlattice that is exactly related to the polarization of 
the full system and does not require choosing an arbitrary 



reference configuration. For each atomic layer, the LP is 
uniquely determined by the spatial locations of ionic and 
WF centers, and it can easily be computed in any first- 
principles code as a post-processing step after standard 
electronic structure calculations. Although this polariza- 
tion is not directly measurable experimentally, we show 
examples in which the LP precisely quantifies polar dis- 
tortions throughout the superlattice, the high resolution 
of the definition being particularly relevant for inspecting 
the behavior of interface layers. Immediate applications 
include modeling of interface effects on total polarization 
of multicomponent superlattices , systematic studies 
of self-poling effects in superlattices j23 , and studies of 
the coupling of phonons to the interfaces. For superlat- 
tices containing magnetic constituents, the spin degener- 
acy assumption can be relaxed, so that the WF centers 
will have additional local spin ordering information. 
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